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Abstract — Excitation of radial oscillations in population I {X = 0.7, Z = 0.02) red 
supergiants is investigated using the solution of the equations of radiation hydrodynamics and 
turbulent convection. The core helium burning stars with masses 8Mq < M < 20Mq and 
effective temperatures T^s < 4000K are shown to be unstable against radial pulsations in the 
fundamental mode. The oscillation periods range between 45 and 1180 days. The pulsational 
instability is due to the K-mechanism in the hydrogen and helium ionization zones. Radial 
pulsations of stars with mass M < 15Mq are strictly periodic with the light amplitude AMboi < 
0™5. The pulsation amplitude increases with increasing stellar mass and for M > 15Mq the 
maximum expansion velocity of outer layers is as high as one third of the escape velocity. The 
mean radii of outer Lagrangean mass zones increase due to nonlinear oscillations by < 30% in 
comparison with the initial equilibrium. The approximate method (with uncertainty of a factor 
of 1.5) to evaluate the mass of the pulsating red supergiant with the known period of radial 
oscillations is proposed. The approximation of the pulsation constant Q as a function of the 
mass-to-radius ratio is given. Masses of seven galactic red supergiants are evaluated using the 
period-mean density relation. 
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1 INTRODUCTION 



Red supergiants are long-period variables with semiregular light variations on a timescale 
> 10^ day. The period-luminosity relation (Glass 1979; Feast et al. 1980) and the linear 
theory of adiabatic oscillations (Stothers 1969, 1972) allow us to suppose that such a type 
of variability is due to radial stellar pulsations. At the same time together with semiregular 
variability some red supergiants exhibit superimposed irregular light variations (Kiss et al. 
2006). The secondary stochastic variability is thought to be due to the large-scale convection 
in the outer subphotospheric layers (Stothers and Leung 1971; Schwarzschild 1975; Stothers 
2010). Red supergiants are also remarkable due to intensive mass loss revealed through a large 
infrared excess indicating dust production in the stellar wind (Verhoelst et al. 2009). 

The period-luminosity relation of radially pulsating red supergiants is used for determina- 
tion of extra-galactic distances and in comparison with Cepheids the red supergiants allow us 
to substantially extend the distance scale due to their higher luminosities (Pierce et al. 2000; 
Jurcevic et al. 2000). Application of the theory of stellar pulsation to the analysis of observed 
variability of red supergiants allows us to verify some conclusions of the stellar evolution the- 
ory in a way similar to that employed earlier for Cepheids. It should be noted also that the 
growing bulk of recent observations indicate that the strong stellar wind of massive late-type 
supergiants is due to nonlinear stellar oscillations (van Loon et al. 2008) 

The nature of radial oscillations in red supergiants is still not completely clear yet. The 
linear analysis of pulsational instability of red supergiants with masses 15Mq < M < 30Mq was 
performed by Li and Gong (1994) and Guo and Li (2002). According to their calculations radial 
oscillations of red supergiants are due to instability of the fundamental mode and, perhaps, the 
first overtone. However the theoretical period-luminosity relation agrees only with fundamental 
mode oscillations. 

Nonlinear radial oscillations of red supergiants were considered only in two studies. In the 
first one (Heger et al. 1997) the authors investigated radial oscillations of the red supergiant 
with mass M = ISM© at the final stage of the core helium burning. In the second work 
(Yoon and Cantiello 2010) the authors investigated pulsational instability of red supergiants 
with masses ISMq < M < 40Mq. It should be noted that in both these studies the self- 
exciting stellar oscillations were treated with modified methods of stellar evolution calculation 
and cfi^ccts of interaction between pulsation motions and turbulent convection were not taken 
into account. 

Below we present results of investigation of nonlinear pulsations of red supergiants obtained 
from the self-consistent solution of the equations of radiation hydrodynamics and turbulent 
convection. The need for such an approach is due to the significant length and mass of the 
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outer convection zone involved in pulsation motions. The treatment of convective heat transport 
uses the solution of the diffusion-type equations for the enthalpy and the mean turbulent 
energy obtained by KuhfuB(1986) for spherically-symmetric gas flows from the Navier-Stokes 
equation. Thus, the results presented below deal with modelling the semiregular variability 
and the secondary stochastic variability is not considered because this problem is beyond the 
approximation of spherical geometry. We consider the stars with masses at the zero-age main 
sequence 8Mq < Mzams < 20Mq and initial fractional mass abundunces of hydrogen and 
elements heavier than helium X = 0.7 and Z = 0.02. 

2 INITIAL CONDITIONS 

The problem of self-exciting stellar oscillations is the Cauchy problem for equations of hy- 
drodynamics with initial conditions corresponding to the hydrostatic and thermal equilibrium. 
In the present study the initial conditions were taken from the evolutionary models of stars at 
the core helium burning. Methods of stellar evolution calculations are discussed in our previous 
papers (Fadeyev 2011a, b). 

The evolutionary tracks in the Hertzsprung-Russel (HR) diagram of population I stars 
{X = 0.7, Z = 0.02) with initial masses Mzams = 10, 15 and 2OM0 are shown in Fig. [1] Final 
points of tracks correspond to helium exhaustion in the stellar core (l^c ~ 10~^). Solid lines 
indicate parts of tracks when the star is the core helium burning red supergiant with effective 
temperature Teg < 4000 K. 

The star leaves the main sequence after core hydrogen exhaustion and crosses the HR 
diagram to lower effective temperatures where becomes the red supergiant with gradually in- 
creasing luminosity. However, from the observational point of view the red supergiants with 
gravitationally contracting helium core are not interesting since this evolutionary stage proceeds 
in the Kelvin-Helmholtz time scale. For example, for the star with initial mass Mzams = lOM© 
the duration of gravitational contraction of the helium core is ~ 10^ years. 

The luminosity ceases to increase when the triple alpha process becomes the main energy 
source in the stellar center. For the star witn initial mass Mzams = lOM© the total duration of 
thermonuclear helium burning is ~ 2.8 ■ 10^ years and in the beginning the star remains the red 
supergiant with luminosity decreasing by a factor of two during ^ 1.2-10® years. The star leaves 
the red supergiant domain when the central helium abundance decreases below Yc ~ 0.52 and 
its evolutionary track loops the HR disgram to effective temperatures as high as T^s ~ lO^K. 
The star becomes again the red supergiant when its central helium abundance decreases to 
Yc ~ 0.05 and the time of helium exhaustion does not exceed 2 ■ 10^ years. Therefore, of 
most interest is the initial stage of helium burning during of which the luminosity of the red 
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supergiant decreases. Evolution of stars with initial masses Mzams < ISM© is nearly the same 
but proceeds in different time scales. 

Evolutionary tracks of the core helium burning stars Mzams ^ ISM© do not loop in the 
HR diagram and all the time remain in the red supergiant domain. For Mzams = SOM© the 
core helium burning proceeds during ^ 8.4 • 10^ years, in the beginning during pa 3.6 • 10^ years 
the stellar luminosity decreasing and then increasing to L ^ IO^Lq. 

All red supergiant evolutionary models used as initial conditions for hydrodynamic compu- 
tations are chemically homogeneous between the inner boundary to the stellar surface. More- 
over,for the given value of Mzams the abundances in the stellar envelope do not change during 
the core hehum burning. This is due to the fact that the size of the outer convective zone is 
maximum at the final stage of gravitational contraction of the helium core just before ignition 
of the triple alpha process. Depending on the mass and luminosity of the star the radius of 
the inner boundary ro ranges within 0.01 < ro/Req < 0.1, where i?eq is the radius of the upper 
boundary of the equilibrium model. 

3 RESULTS OF HYDRODYNAMIC COMPUTATIONS 

The method for the self-consistent solution of the equations of radiation hydrodynamics and 
turbulent convection is described in our previous paper (Fadeyev, 2011b), so that below we only 
discuss the results obtained. Computations were carried out with the number of Lagrangean 
mass zones SOO < TV < 10^. To be confident that the solution is independent of the inner 
boundary radius ro and the number of Lagrangean zones N some hydrodynamic models were 
computed with several different values of these quantities. 

Our hydrodynamic computations show that red supergiants with initial masses 8Mq < 
Mzams < 2OM0 are unstable against radial oscillations in the fundamental mode. However 
depending on the value of Mzams hydrodynamic models demonstrate different behaviour both 
during the growth of instability and after the limit amplitude attainment. In particular, in stars 
with initial mass Mzams < ISM© the oscillation amplitude is always enough small, whereas for 
Mzams > ISM© nonlinear effects play a substantial role. 

The main properties of hydrodynamical models are summarized in Table 1, where the 
bolometric luminosity L and the effective temperature Tcs correspond to the initial equilibrium 
model, Yc is the fractional helium abundance in the stellar center, 11 and Q are the pulsation 
period and the pulsation constant in days, rj = UdlnEK/dt is the growth rate of the pulsation 
kinetic energy £Jk- Reciprocal of this quantity equals the number of pulsation periods during 

which the kinetic energy increases by a factor of e = 2.718 Ratios of the outer boundary 

mean radius (R) to the equilibrium radius of the model Req are given in the last column of 
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Table 1 and show the role of nonlinear effects after the limit amplitude attainment. 

The oscillation amplitude of red supergiants Mzams = lOM© at the top of the evolutionary 
track is enough small and the relative radial displacement of the upper boundary is AR/Rcq = 
0.06. The oscillation amplitude gradually increases with decreasing luminosity and in the point 
with minimum luminosity before the loop in the HR diagram AR/Req = 0.13. Enhancement 
of the radial oscillation amplitude with decreasung luminosity is illustrated in Fig. [2] where 
variations of the velocity at the upper boundary U and the bolometric magnitude Mboi are 
shown for L = 2 ■ 10^ 1.5 ■ 10^ and lO^L©. 

Fairly good repetition of small amplitude oscillations allows us to calculate the mechanical 
work done by Lagrangean mass zones and thereby to evaluate their contribution into excitation 
or damping of instability. The radial dependence of the mechanical work f PdV, where V is the 
specific volume and P is the sum of gas, radiation and turbulent pressure, is shown in Fig. [3] 
for the hydrodynamical model Mzams = IOMq, L = 1.5 ■ IO^Lq. The region of instability 
excitation PdV > 0) encompasses the layers with temperature 1.2 ■ 10^ ^ T < 4 ■ lO'^K 
corresponding to the hydrogen and helium ionization zones. In deeper layers (T > 4 ■ lO'^K) 
with fully ionized helium the pulsational instability is damped PdV < 0). 

To understand the physical mechanism of excitation of pulsational instability let us consider 
variations of the gas density p, temperature T, opacity k and luminosity Lrad = 47rr^Fi.ad, where 
Frad is radiative flux, in Lagrangean mass zones of the hydrodynamical model. In Fig. Hl^a) we 
give the plots of relative variations Sp/p, 6T/T and 6k/ k in the layers of fully ionized helium 
with temperature ranging within 4.9 ■ 10^ K < T < 5.6 ■ lO'* K. For the sake of convenience, the 
plots are arbitrarily shifted along the vertical axis. Coincidence of the maxima of density and 
temperature variations indicates that oscillations are nearly adiabatic. Decrease of opacity at 
maximum compression damps the instabihty because, as seen in Fig. Hl^b), heat losses due to 
radiation reach their maximum. 

Variations of same quantities for the layer with temperature 1.4 ■ 10^ K < T < 1.7 ■ 10^ K 
corresponding to partial helium ionization are plotted in Fig.O Substantial phase shifts between 
maxima of density and temperature indicate large nonadiabaticity of pulsation motions, whereas 
the delay of the maximum temperature with respect to maximum compression is the cause of 
the positive mechanical work. Absorption of heat during compression is due to increase of 
opacity and it is accompanied by decrease of the radiative flux. 

Thus, damping of oscillations in the layers of fully ionized helium and driving of pulsational 
instability in the hydrogen and helium ionization zones are due to the K-mechanism, because 
effects of heat gains and losses are connected with absorption and emission of radiation. Low 
rates of the instability growth (77 ~ 10^^) and the small limit cycle amplitude are due to the 
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small fraction of radiation in the total energy flux (Lj-ad ^ 10~^Lr). Driving of pulsational in- 
stability at so small radiation fluxes is due to the large amplitude of total luminosity variations. 
As seen in Fig. [6l the amplitude of luminosity variations is largest in vicinity of the helium 
ionization zone. 

Pulsational instability of red supergiants increases with increasing initial stellar mass and 
for MzAMS > 15Mq nonlinear effects become significant. In Fig. [7] we give the temporal 
dependences of the kinetic energy Ek of the pulsating envelope and the radius of the upper 
boundary R in units of the equilibrium radius i?eq for the red supergiant model with initial mass 
A^zAMS = I6M0 and luminosity L = 8.5 ■ lO^L©. Compared to less massive supergiants this 
model demonstrates the growth rate which is higher by an order of magnitude, whereas after 
the attainment of the limiting amplitude the mean radii of outer Lagrangean mass zones exceed 
their equilibrium values by nearly one third. The amplitude ceases to grow in a transitional 
process encompassing roughly a dozen of oscillation periods (20 < t/U < 30). During this time 
interval the sources of instability in the hydrogen and helium ionization zones are balanced by 
damping of instability in layers of fully ionized helium and shock radiative losses in the layers 
above the photosphere. 

4 PERIOD-LUMINOSITY RELATION 

In the HR diagram red supergiants occupy the domain with relatively narrow effective 
temperature range (3000 K < T^s < 4000 K), so that similar to Cepheids they exhibit cor- 
relation between the equilibrium luminosity L and the period of radial pulsations 11. From 
observations such a correlation is known as the period-luminosity relation. The theoretical 
period-luminosity relation obtained in the present study is shown in Fig. |H] for models of three 
evolutionary sequences with initial masses Mzams = 10, 15 and 20Mq. Hydrodynamical mod- 
els of red giants at the evolutionary stage of decreasing luminosity are shown by filled circles 
and the dashed lines show the linear fits for each evolutionary sequence. 

For stars with initial mass Mzams = IOMq the evolutionary track loops the HR diagram 
during helium burning and three models shown in Fig. [8] by triangles correspond to the initial 
part of the loop with effective temperatures T^s < 4000K. Therefore, one of the causes of the 
dispersion of points in the empirical period-luminosity diagram is that some stars with masses 
M < I5M0 leave the red supergiant domain and others return to it. 

The chemical composition of outer layers involved in pulsation motions does not change, 
whereas effects of mass loss for M < I5M0 are insignificant. Therefore, red supergiants with 
almost exhausted helium in the stellar center obey the same period-luminosity relation as stars 
in earlier phases of helium burning. This is illustrated in Fig. [8] where two models of stars 
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^ZAMS = IOMq with central helium abundances Y^. = 2.4 ■ 10"'^ and 1.1 ■ 10~^ are shown by 
open circles. 

As can be seen in Fig. [HI the dispersion of the common correlation between the period 
and luminosity of red supergiants is due to dependence of the both luminosity and period on 
the stellar mass. Therefore, one of the causes of dispersion on the empirical period-luminosity 
diagram is due to different masses of observed stars. It should also be noted that the mass- 
luminosity relation of red supergiants and, therefore, the period-luminosity relation, depend on 
convective overshooting and mass loss during the preceeding evolution. An important role in 
the both mass-luminosity and period-luminosity relations belongs also to abundances of heavy 
elements Z. These effects, however, were beyond the scope of the present study. 



The equlibrium luminosity of the red supergiant during helium burning changes within the 
ranges that depend on the initial stellar mass. For example, in stars with Mzams ~ lOM© the 
luminosity decreases by a half, whereas in red supergiants with initial mass Mzams = SOM© the 
maximum to minimum luminosity ratio decreases to ~ 1.6. The period of radial pulsations 11 
changes simultaneosly with equilibrium luminosity L. Evolution of red supergiants Mzams < 
I2M0 between the upper and lower luminosity limits is accompanied by the change of the 
pulsation period by a factor of three. The maximum to minimum period ratio decreases to a 
factor of two for Mzams = 2OM0 . 

This property of red supergiants is illustarted in Fig. |9] where for stars with initial masses 
8Mq < Mzams < 2OM0 we show the period-mass diagram. Hydrodynamical models of stars 
at the top of the evolutionary track are shown by filled circles. Open circles indicate the red 
supergiant models with lower luminosity. The diagram in Fig. |9] takes into account effects of 
mass loss and along the vertical axis we give the mass values M of evolving stars. Evolution of 
the red supergiant corresponds to the displacement on the diagram from right to left and then 
in the opposite direction. For models Mzams = IOMq and 2OM0 this displacement is shown 
by arrows. 

The period-mass diagram in Fig. [9] demonstrates the existence of the limited area of mass 
and period values. The borders of allowed masses and periods of radially pulsating red super- 
giants can be approximately fitted as 



and in Fig. IH] they are shown by dashed lines. Thus, from the observational estimate of the 
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period 11 relations ([T]) allow us to evaluate the upper and lower mass limits of the red supergiant. 
For periods 11 < 300 the uncertainty of such an estimate is about a factor of ~ 1.5. 

6 PULSATION CONSTANT 

The more exact value of the red supergiant mass can be obtained from the period-mean 
density relation 

n = Q {R/Ref (M/Mq)-'/^ (2) 

provided that the pulsation period 11 and the mean stellar radius R are known from observations. 
The pulsation constant Q is obtained from the theory of stellar pulsation and in some cases 
can be expressed as a function of the stellar mass M and stellar radius R. Substitution of this 
expression into the period-mean density relation ([5]) allows us to obtain the explicit expression 
for the mass of the pulsating star. 

The pulsation constants obtained in our hydrodynamical calculations of red supergiants 
with initial masses 8Mq < Mzams < 20Mq and pulsation periods 45 < 11 < 1180 are shown 
in Fig. [TD] as a function of mass-to-radius ratio / = {M/Mq)/{R/Rq). The linear fit of the 
pulsation constant is written as 

log g = -2.288 -0.778 log/ (3) 

and is shown in Fig. [10] by the dashed line. 

Masses M of seven galactic red supergiants evaluated from substitution of ([3]) into the 
period-mean density relation are given in Table 2. The periods 11 are taken from the 
General Catalogue of Variable Stars (Samus et al. 2011). The mean stellar radii were evaluated 
by Levesque et al. (2005) and Josselin and Plez (2007). In last two columns of Table 2 we give 
the lower Ma and upper mass limits derived from ([1]). 

Unfortunately, the existing estimates of mean radii of red supergiants are still highly un- 
certain. For example, the uncertainty in the effective temperature is ~ 25% (Josselin and Plez 
2007), so that the uncertainty in the mean radius is as high as ~ 60%. Therefore, the case 
when the stellar mass determined from the period-mean density relation is beyond the interval 
[Ma, Mh] should not be considered as a contradiction. For example, if we adopt that the radius 
of AD Per is larger by 20% (R = 548/?©) then the stellar mass is M = 12.3M0, that is within 
ranges given by 
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7 CONCLUSION 



Given in the previous section estimates of masses of seven galactic red supergiants allow us 
to conclude that the theory of stellar evolution is in an agreement with observational estimates 
of stellar radii. To compare more stars with the theoretical period-luminosity relation one 
should consider pulsational instability of red supergiants in the wider interval of initial masses 

-^ZAMS- 

For more detailed theoretical period-luminosity relation one should consider the role of some 
parameters used in evolutionary computations. One of them is the overshooting parameter. In 
the present study the evolutionary computations were done for the ratio of the overshooting 
distance to the pressure scale height low/ Hp — 0.15. The need to know the role of this parameter 
is due to the dependence of the mass-luminosity relation of helium burning stars on convective 
overshooting. 

In stars with masses M > 20Mq effects of mass loss during the red supergiant evolutionary 
stage become significant. In the present study the evolutionary calculations were done with 
mass loss rates by Nieuwenhuijzen and de Jager (1990) however determination of the mass 
loss rate M as a function of fundamental stellar parameters remains disputable (Mauron and 
Josselin, 2011) Therefore, one should employ parametrization of the expression for M and 
consider the mass-luminosity and period-luminosity relations as a function of this parameter. 

Another parameter which significantly affects the period-luminosity relation of red super- 
giants is the mass fraction abundance of heavy elements Z. Of special interest is the period- 
luminosity relation for Z — 0.008 which is typical for the Large and Small Magellanic Clouds. 
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Table 1: Hydrodynamical models of red supergaints. 



Mzams/Mq 


M/Mq 








n, 


Q, 


V 


{R)/Req 






10^ 


K 




day 


day 






20 


19.53 


14.029 


3320 


0.975 


1178 


0.1365 








18.36 


9.380 


3478 


0.746 


608 


0.1058 


0.27 


1.24 




17.47 


8.782 


3522 


0.480 


555 


0.1028 


0.30 


1.11 


15 


14.80 


7.176 


3370 


0.968 


651 


0.1131 


0.18 


1.26 




14.12 


4.454 


3556 


0.680 


310 


0.0883 


0.10 


1.12 




13.69 


4.074 


3620 


0.379 


262 


0.0830 


0.11 


1.11 


10 


9.88 


2.055 


3570 


0.970 


203 


0.0867 


0.01 


1.15 




9.76 


1.507 


3672 


0.816 


127 


0.0747 


0.03 


1.10 




9.66 


1.102 


3870 


0.523 


67 


0.0615 


0.03 


1.04 


9 


8.90 


1.452 


3610 


0.973 


154 


0.0849 


0.02 


1.16 




8.76 


0.668 


3958 


0.580 


45 


0.0574 


0.02 


1.02 


8 


7.02 


0.962 


3690 


0.972 


102 


0.0765 


0.02 


1.13 
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Table 2: Masses of galactic red supergiants. 





n, day 


R/Rq 


M/Mq 


MJMq 




SU Per 


533 


780 


17.1 


14.1 


17.1 


W Per 


485 


620 


12.2 


13.6 


16.6 


V602 Car 


635 


860 


17.7 


15.0 


17.9 


AD Per 


362.5 


457 


8.9 


12.2 


15.4 


FZ Per 


184 


324 


8.2 


9.5 


12.8 


RWCyg 


550 


676 


12.9 


14.2 


17.2 


SW Cep 


70 


234 


9.8 


6.7 


9.8 
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Figure captions 

Fig. 1. Evolutionary tracls of Population I stars with initial masses Mzams = 10, 15 and 20Mq 
in the HR diagram. Parts of tracks corresponding to the stage of the red supergiant are 
shown in solid lines. 

Fig. 2. Variations of the gas flow velocity at the upper boundary U (a) and bolometric magnitude 
M]joi (b) in red supergiants with initial mass Mzams = lOM© and luminosity L = 2 ■ 10^ 
(solid lines), 1.5 ■ lO'^L© (dashed lines) and lO^L© (dot-dashed lines). 

Fig. 3. The radial dependence of the mechanical work done by a Lagrangean mass zone over the 
pulsation period 11. 

Fig. 4. - Variations of gas density p (solid line), temperature T (dashed line) and opacity k, 
(dotted line) in the layer of fully ionized helium; (b) - variations of radiative luminosity 
in units of the total equilibrium luminosity Lq. 

Fig. 5. Same as Fig. IHbut for the layer with partial helium ionization. 

Fig. 6. The amplitude of variations of the total luminosity expressed in units of the equilibrium 
luminosity Lq versus the radial distance from the stellar center. 

Fig. 7. Kinetic energy of pulsation motions -Ek (a) and the upper boundary radius R (b) as a 
function of time t for the hydrodynamical model Mzams = I6M0, L = 8.5 ■ IO'^Lq. 

Fig. 8. The period-luminosity diagram for red supergiants with initial masses Mzams = 10, 15 
and 2OM0. Hydrodynamical models are represented by filled cicrles. In triangles are 
shown hydrodynamical models of stars that leave the red supergiant domain. In open 
circles are represented the hydrodynamical models of stars with central helium abundance 
Fc < 2.4 ■ 10-3. 

Fig. 9. The period-mass diagram of red supergiants 8M0 < Mzams < 2OM0. Hydrodynamical 
models of stars with maximum and minimum equilibrium luminosity are shown by filled 
and open circles, respectively. Arrows indicate the direction of evolutionary change of the 
pulsation period H of stars Mzams = IOM0 and 2OM0. The region of allowed values of 
the radial pulsation periods is limited by dashed lines. 

Fig. 10. The pulsation constant Q of red supergiants 8M0 < Mzams < 2OM0 as a function of 
mass-to-radius ratio / = (M/Mq)/ {R/ Rq). Hydrodynamical models of red supergiants 
are shown by filled circles. 
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4.4 4.0 3.6 

log T,,f 

Figure 1: Evolutionary tracls of Population I stars with initial masses Mzams = 10, 15 and 
2OM0 in the HR diagram. Parts of tracks corresponding to the stage of the red supergiant are 
shown in solid lines. 
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Figure 2: Variations of the gas flow velocity at the upper boundary U (a) and bolometric 
magnitude Mboi (b) in red supergiants with initial mass Mzams = lOM© and luminosity L — 
2 • 10^ (solid hues), 1.5 • IO^Lq (dashed lines) and IO^Lq (dot-dashed hues). 
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Figure 3: The radial dependence of the mechanical work done by a Lagrangean mass zone over 
the pulsation period 11. 
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Figure 4: () - Variations of gas density p (solid line), temperature T (dashed line) and opacity 
K (dotted line) in the layer of fully ionized helium; (b) - variations of radiative luminosity in 
units of the total equilibrium luminosity Lq. 
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Figure 6: The amplitude of variations of the total luminosity expressed in units of the 
equilibrium luminosity Lq versus the radial distance from the stellar center. 
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Figure 8: The period-luminosity diagram for red supergiants with initial masses AfzAMS = 10) 
15 and 20Mq. Hydrodynamical models are represented by filled cicrles. In triangles are shown 
hydrodynamical models of stars that leave the red supergiant domain. In open circles are 
represented the hydrodynamical models of stars with central helium abundance Yc < 2.4 • 10~^. 
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Figure 9: The period-mass diagram of red supergiants SM© < Mzams < 2OM0. Hydrody- 
namical models of stars with maximum and minimum equihbrium luminosity are shown by 
filled and open circles, respectively. Arrows indicate the direction of evolutionary change of the 
pulsation period IT of stars Mzams = IOMq and 2QMq. The region of allowed values of the 
radial pulsation periods is limited by dashed lines. 
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Figure 10: The pulsation constant Q of red supergiants 8Mq < Mzams < 20Mq as a function 
of mass-to-radius ratio / = [M / Mq) / {R/ Rq). Hydro dynamical models of red supergiants are 
shown by filled circles. 
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